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Atomic-scale logic and the minimization of heating (dissipation) are both very high on the agenda for 
future computation hardware. An approach to achieve these would be to replace networks of transistors 
directly by classical reversible logic gates built from the coherent dynamics of a few interacting atoms. 
As superpositions are unnecessary before and after each such gate (inputs and outputs are bits), the de¬ 
phasing time only needs to exceed a single gate operation time, while fault tolerance should be achieved 
with low overhead, by classical coding. Such gates could thus be a spin-off of quantum technology much 
before full-scale quantum computation. Thus motivated, we propose methods to realize the 3-bit Toffoli 
and Fredkin gates universal for classical reversible logic using a single time-independent 3-qubit Hamil¬ 
tonian with realistic nearest neighbour two-body interactions. We also exemplify how these gates can be 
composed to make a larger circuit. We show that trapped ions may soon be scalable simulators for such 
architectures, and investigate the prospects with dopants in silicon. 


Introduction:- Power dissipation has become a seri¬ 
ous obstacle to packing more transistors per unit area so 
that the exponential rise of computational power with time 
(Moore’s law) may be continued. Heating in a chip is pro¬ 
jected to reach 200Wcm“^ by 2020 ll]] by the International 
Technology Roadmap for Semiconductors (ITRS). The 
search for less dissipative alternatives is on with sugges¬ 
tions such as molecular electronics ||2], spin-wave compu¬ 
tation Q, u iag netic and quantum dot cellular automata 0], 
DNA logicand superconducting logic in cryogenic tem¬ 
peratures 0, which are by no means exhaustive. 

Aside from minimizing dissipation, another driving fac¬ 
tor for contemporary computer technology is the atomic 
scale storage of information illTii and atomic scale 
logic J^. However, to our understanding, they do not yet 
aim to exploit the dynamics of highly isolated systems for 
computation in the same sense as the “friction free” bil¬ 
liard ball computer of Fredkin and Toffoli ifT^ . Such dy¬ 
namics is well approximated in the systems being devel¬ 
oped for quantum technologies where quantum coherence 
is preserved by a high isolation. In fact, the energy dissi¬ 
pation time-scale Ti can be exceptionally high, and even 
the dephasing time T 2 is fairly high. It is thereby worth 
studying whether the huge development towards quantum 
computation can, on the way to that grand goal, also pro¬ 
vide a minimally dissipative, as well as miniaturized, hard¬ 
ware for classical logic. Any reliability sacrificed by going 
to atomic bits is not particularly new at this scale, and, is 
present even at the nano-scale, and ingenious ways of using 
low-reliability devices is topical ifTl- 131. 

Motivated by the above, here we investigate whether an 
unmodulated minimal widget of 3 permanently interacting 
qubits, with each qubit encoding a microscopic bit, can act 
as a logic gate for classical computation. The aim is to (a) 
use coherent dynamics (to avoid dissipation and heating), 
(b) use permanent nearest neighbour two-body couplings 


of similar strength (to keep things realistic for a structure of 
proximal spin qubits, for example), (c) use a single “time- 
independent” Hamiltonian to accomplish the entire gate 
and finally, (d) avoid any auxiliary systems, hybrid systems 
or additional levels aside those of the relevant qubits. In 
practical implementations, additional regular pulsings may 
however be required for reducing decoherence through dy¬ 
namical decoupling. We suggest placing these gates next 
to each other spatially to compose a classical circuit (we 
also exemplify this composability). The whole classical 
circuit will then simply be a 2D pattern of the widgets im¬ 
plementing the fundamental gates, where each such gate is 
implemented by a quantum evolution. However, the quan¬ 
tum state is allowed to decohere “before” and “after” each 
gate. This is acceptable for a classical circuit as the inputs 
and outputs of each gate are classical bits, and no superpo¬ 
sition has to be maintained between the end of one basic 
gate and the start of another. Unlike the case for quantum 
computation, in the space of time between the gates, the 
correction of dephasing errors are not necessry. The use of 
classical codes therefore suffices to keep track of the errors 
within the computation. 

As coherent dynamics is not only non-dissipative, but 
also reversible, we will aim to build 3-bit Toffoli and 
Fredkin gates, which enable reversible classical compu¬ 
tation 113 B. Reversible log ic avoids heating due to 
the erasure of information 111^ . Quantum computation 
and quantum error correction has already motivated the 
implementation of Toffoli gates Although these 

gates can be decomposed into 2-qubit gates and local uni- 
taries id, such implementations require at least five 2 
qubit gates lE5l . To achieve simplifications beyond this 
remit (including single pulse or “single-shot” implementa¬ 
tions), nearly all schemes and actual implementations have 
variously used auxiliary modes aside the relevant qubits, 
such as the cavity mode in hybrid qubit-resonator systems 
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EEMl or the motional modes in ion traps 
or auxiliary levels outside the space of quhits 1251 
Where solely quhits have been used, either multiple pulses 
or non-uniform and long-range couplings (departing 
from aim (h)) are generically present as in the NMR and 
other literature j^BlIl . Crucially, the key question of the 
“possibility” and an analytic expression of the “accuracy” 
of gates in the simplest setting: 3 qubits and a time inde¬ 
pendent realistic nearest neighbour Hamiltonian, remains 
open. In fact, the conjuction (a)-(d) should be impossible 
as the generation of the unitary operations corresponding 
to the ideal Fredkin and Toffoli gates with a single (time- 
constant) 3 qubit Hamiltonian seems to necessitate unreal¬ 
istic 3-body interactions |l32, 33] cr^crj ■ cTk and 
respectively. However, it has gone unnoticed that when we 
lower our aim from quantum to classical computation, i.e., 
when the relative phases between the computational basis 
states are irrelevant, and when approximate gates with low 
errors could be useful, then classical Fredkin and Toffoli 
gates of useful accuracy become feasible with 3 qubit real¬ 
istic nearest neighbour time-independent Hamiltonians. 

Toffoli gate:- A classical Toffoli gate flips the target bit 
when both the control bits are in the logical state 1 (a 
quantum Toffoli is a unitary operator that additionally en¬ 
codes a specific phase relationship between distinct quan¬ 
tum states). We start by describing how a structure of 3 
permanently Ising coupled spins can be used to implement 
an approximate classical Toffoli gate and investigate how 
good the approximation can be. The gate is switched on 
by applying a transverse held to the target qubit (qubit 2 in 
this case, see fig. [T]). With this held switched on, the gate 
is performed through the time-independent Hamiltonian 
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where J, 


ZZl^j 


and H are in frequency units (unless oth¬ 
erwise specified we will use fhese units). For the Toffoli 
gate we require -^|l)(|0)±e*‘^|l))|l)tobe eigenstates of 
FfxoF- By applying F^tof on states of this form, it can be 
confirmed fhat 0:2 = 2Jzz achieves the desired eigenstate 
with ^ = 0, while tci and can remain arbitrary. The 
eigenstates and energies are then 
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FIG. 1. Setup for creating a 3-qubit Toffoli gate. 


where i,j € {0,1}, du = 0, dot = <^10 = W2/O, doo = 
2a;2/n, and are normalising factors. The |101) o 
|111) swapping will occur when exp [—fiFTori] |101) = 
6*^1 111), which occurs at a time t = Tn = (2n -f 
l)7rh/\Ey, - Eyy = {2n -f l)7r/0 (where n is an inte¬ 
ger, and assuming without loss of generality that Q is pos¬ 
itive). In general, the evolution of the arbitrary computa¬ 
tional basis states in this time is captured by the fideli¬ 
ties fimn^xyz ■= |/mn). As [FFtof,^-^] = 

[FFxofi O'!] = 0 only the following fidelities are relevant 
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where a,b,c G {0,1}, b := F © 1, and 4>ac = 
( 2 n+l).((-ir^, + (-l)-a; 3 ) _ 1/101^111 1 = 1 by 

our choice of . To realise a Toffoli gate we further require 
that l/ioj ^ fiijl = 0 for \ij) y |11) so that for these fi¬ 
delities, the phase inside the sine function in (|4| must be an 
integer multiple of tt. This leads us to 
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where mi, m2 are non-zero integers, which in turn implies 
16mi — 4m2 = 3(2n+1), where the left hand side is even, 
and the right side is odd. Thus no choice of a;2, O gives a 
perfect Toffoli (a price to pay for the simplicity of FFtof)- 
However, we can And parameters which achieve an ap¬ 
proximate gate. Assuming n = 0 for the shortest pos¬ 
sible gate time and further assuming that the first sube¬ 
quation of Eq.® is exact with large mi one finds m2 ~ 
2mi from the second subequation of Eq.(|5]l. With this 
choice of parameters, the phases inside the sine func¬ 
tions in dH are all either multiples of tt or approximately 
so up to order 1/mi. To evaluate how close this ap¬ 
proximate Toffoli is to the exact Toffoli, we use the pro¬ 
cess trace distance, which for a 3-qubit system is defined 
as [|34j] T)pno tr\x{U) - x{T)\ /16 where x{M) mn 
tr(Aj„M)(tr(At M))*, and ^ complete or¬ 

thogonal basis which satisfies tr(AT = 5nm, T is an 
ideal Toffoli gate, U is the gate we can achieve with the 
above setup, and 1X1 indicates the matrix norm. As we 
are interested in creating a “classical” gate, we ignore any 
phases and define such that {y\U\f\\x) = \fx^y\ so 
this will only measure how close we are to a Toffoli gate 
apart from local operations. gives an upper bound on 
the average probability pe that the gate fails ydj], so 
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In summary, we can achieve an approximate classical Tof- 
foli gate with average failure error of Eq.®, if = ^ 2/2 
and a; 2 /n = — 1 , with mi large. 

Fredkin Gate:- We now consider creating a classical 
Fredkin gate (controlled-SWAP), using quantum Ising and 
Heisenberg interactions. We consider the Hamiltonian 

f^FRED = "^^2 ■ cr 3 H—^CTicrl + —-—, (7) 

3 

where quhit 1 is the control quhit, and quhits 2 and 3 are 
to he swapped (see fig. |2l). The intuition is that the swap¬ 
ping induced hy the Heisenberg interaction between qubits 
2 and 3 will only occur when qubit 1 is in a state that makes 
the energy splitting of qubits 2 and 3 match. For swap be¬ 
tween qubits 2 and 3, we need to choose parameters such 
that states of the form -^(IllO) ± | 101 )) are eigenstates of 
^FRED> which requires J^z = u ;2 — ^ 3 . With these parame¬ 
ters, {| 100 ), | 011 ), | 000 )} are all eigenstates, and the 

only eigenstates of the Hamiltonian which are not compu¬ 
tational basis states are 


J 12 = J 23 = Jzz = 27r X 9.98kHz, J 13 = 27r x 7.07kHz. 
The extra J 13 coupling introduces extra phases to the gate 
but this has no effect if the gate is used for classical compu¬ 
tation. The (jf field is achieved by applying a near-resonant 
microwave pulse leading to the trapped ion Hamiltonian 
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where a;9 is the qubit energy splitting of ion j. 
can transform 77^*^ to a frame rotating with the operator 
~ ^^ 2 )’ where S is the detuning of the mi¬ 
crowave field from the resonant frequency of ion 2. States 
in this rotating frame evolve according to 


We 


= 


J: 


9 / ^ ^3 ^3 + 1 


+ cos{u3R)al 


J = 1 


Choosing 5 = 2 Jzz and iOx = UJ 2 — S, and applying iden¬ 
tity [13] of the supplementary material, we have 


= (| 110 ) ± \Wl))/V2, = {^ 1010 ) + (^^2 

-CU3 ± v'(w 2-CU3)2 + j2)|001)}/Af5to (8) 


where is a normalising factor. The eigenener- 

gies of these states are = — (wi + J =F 2 J) /2 

and E^iq/H = — J ± 2^Jzz + /2 respectively. 

The swap |110) o |101) is complete at a time = 
(2n -I- l) 7 rfi/|i 7 ^jQ — + l) 7 r/ 2 J (with n G 

{0,1,2,...}, and assuming J > 0, without loss of general¬ 
ity). The fidelity for swapping |010) fG |001) at time 
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be zero, we need J^ (m^/( 2 n -|- 1 )^ — l) = [u32 — 


where m is an integer greater than 1 , and 
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Possible realizations:- In implementations, the gate op¬ 
eration times, ^ 1/Q for the Toffoli, and ^ 1/J for the 
Fredkin, have to be smaller than the dephasing times. We 
discuss two possible realizations. 

Trapped Ions:- To engineer the Hamiltonian iTroF we ex¬ 
ploit the fact that an axial magnetic field gradient realises 
Ising couplings ll^ . Consider three ^^^Yb“'“ ions in a lin¬ 
ear Paul trap with secular frequency a = 2Trx lOOkHz. The 
qubits are encoded in the ^S' 1/2 l'{) = 1 ^ = 0 ,mi? = 0 ) 
and It) = \F = 1, mp = 1) states which are separated by 
approximately 12.6GHz. For dzBj = 250Tm“^, one has 
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FIG. 2. Setup for creating a 3-qubit Fredkin gate, using Ising and 
Heisenberg coupling. 
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Thus with Q = 27r X 1.8kHz and 


n 


2(^2 ‘^Jza 

12.6GHz, 


w TTtof. while the systematic gate error from Eq.® 


is ^ 0.02. The gate time is ^ 0.3ms, giving an error due 
to decoherence of around (1 — ^ 0.03 with 

T 2 ~ 10ms ll^ using dynamical decoupling). In addi¬ 
tion there may be additional heating due to proximity of 
the ions to the electrode surface, which we estimate from 
the results in ll37ll . For a cryogenically cooled surface trap 
with an ion-electrode distance of 160 pm we estimate an 
additional decoherence rate of 50 Hz, which results in an 
error of around 0.01. Overall we therefore expect an aver¬ 
age gate with error ^ 0.04, which should be good enough 
for error-free classical circuitry ll38l - l4lll . 

Bismuth Donors in Silicon:- We propose placing the 
donor atoms ll42l - l^ close to each other so that their elec¬ 
tronic spins are permanently coupled by isotropic Heisen¬ 
berg interactions. To engineer Hthp, the nuclear spins are 
prepared in different states {/f, /f, If }, resulting 

in different hyperfine couplings (with the very high cou¬ 
pling strength A = 1.475 GHz of Bi) at each site. Starting 
from a nuclear spin polarised sample S we can flip the 
nuclear spin from | to — | in 9 steps of ^10/rs each. The 
nuclear spins are stable for hours so this process 

could be done once before many operations of the gate. We 
use a magnetic field utp » A to ensure the nuclear spin 
does not evolve. The Hamiltonian of the donors is then ll53ll 
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where S'" = ^cr" with a = x,y,z are the Pauli matrices 
for electronic spin. An AC field of strength Q is applied 
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on qubit 2 in the x-direction (which could be applied glob¬ 
ally since the hyperfine splittings are different) so that the 
Hamiltonian + Qcos{uJxt)S 2 acts on the 

donors. Setting J 12 = J 23 = J,^x = wl + “ J 

and transforming to a frame rotating with the opera¬ 
tor ujl ELi Jsi + Yfn=i ^InSn results in (using 

identities [TdIfTdl of the supplementary material) 

_ ^^Ovi^Jzz = J,t 02 = 2J,0Jl = LO3 = 0) 

^ “ 2 

/A j n \ 

^ [^2 ^ 2(a;L - 2J + 2l /|)) ' 

Setting II = |, /| = — |, and /| = |, and with 

0 = IMHz, J = 30MHz, both the error term in and 
the systematic gate error term in Eq.® are ^ 10“^. The 
gate time is 2ps (thereby allowing the bandwidth of the AC 
pulse to significantly exceed the ^ 2kHz linewidths seen in 
experiments fs^l) so that the errors due to decoherence are 
roughly 1 — ^ 10“® (as T2 ~700ms in isotopi- 

cally pure silicon jSSh ). 

The Fredkin gate can be implemented solely with . 
Transforming to a rotating frame with the operator 

ujl En=i Sn+AIISl + Al2{S2+SI) and using the iden¬ 
tity [T 4 ]of the supplementary material gives 

~ = Jl2,J = ^23,^2 

- ujs = 2A{II-II)) + 0{Ju/\AI!-AII\). (11) 

Setting /f = — |, /2 = |i -fs = I, Fredkin gate conditions 

J12 = 2(A/| - All) = J23 ^j^^-T can be met 

by J23 = ^12(1 + 10“®) (n = 675 , m = 2340 ) with a 
resulting gate time ^ 0 . 5 /rs = 10“^r2 so that decoherence 
is negligible. The error term in Eq. lfTTI) is ~ 0.22, but could 
be minimized further to ^ 0.07 by techniques mentioned 
in the supplementary material. 

Composability:- To exemplify circuit building, we show 
how a half-adder (fig. [3] a)) 115^ can be implemenfed wifh 
fhe arrangemenf of Ising coupled spins shown in fig.[3]b). 
The Toffoli gale wifh qubil 2 as largel is implemented using 
TfxoF- Then we wanl lo apply fhe conlrolled-NOT irrespec- 
live of fhe slate of qubil 2. Two successive pulses on qubil 3 
of frequencies a;+ and cc-, where uj± = ws — J13 ± J23 im- 
plemenl Iwo successive condilional flips of qubil 3 accord¬ 
ing lo when qubils 1 and 2 are in a | 10 )x 2 and |ll)x 2 re- 
speclively. Each of Ihese pulses implemenl differenl lime- 
independenl Hamiltonians in appropriate rolaling frames. 
Thus, 3 successive lime independenl Hamiltonians imple- 
menl a half adder. In general, a Toffoli gale can be applied 
on a sel of qubils (say, 1, 2 and 3 wilh 2 as largel) of Ihe 
Ising coupled array depicted in fig.[3]c), by pulses of 4 fre¬ 
quencies to flip qubil 2 irrespeclive of Ihe slate of Ihose 
neighbours (say, A and B) thal do nol lake part in Ihe gale. 
White our pulsing is similar to tools in liquid slate NMR 


il was nol apparenl to date that the gates possible 
in the simplest of settings are approximate, inequivalent to 
the unitary operations corresponding to Fredkin and Toffoli 
gates (to achieve those further “non-local” gates are neces¬ 
sary), and do not require long-range couplings. 



FIG. 3. a) A half-adder circuit, b) Setup for creating a half adder, 
using two pulses and with J 12 = J 23 A Ji3- c) Using selec¬ 
tive addressing on arrays of qubits, general computations can be 
achieved. 

Reliability:- Our physical realizations have errors (as in 
any nano-scale logic, including scaled CMOS). However, 
reliable classical computation with faulty components is 
possible with a constant overhead llH 40|l a s our error 
rates are below the 1/6 required threshold 11^ . One can 
use measurements between gates and a classical resetting 
of bits (e.g., parity protected gates llssi] ). Besides, applica¬ 
tions such as image processing tolerate more noise lll3l] . If 
we automate error correction with faulty gates such as in 
“quantum” error correction, thresholds 10 ^ — 10 ^ are 
however obtained (see Supplementary Material J, which 
are still met by the donor based implementation of Hjqf- 

Conclusions:- We have demonstrated how the classi¬ 
cal Toffoli and Fredkin gates can be achieved by realis¬ 
tic 3 qubit time-independent Hamiltonians. This focus 
on simplicity (e.g., time-independence, no auxiliary sys- 
tems/levels) as opposed to fidelity lEsll . stems from aim¬ 
ing to build low dissipation atomic-scale classical logic. 
Targeting classical gates helps us circumvent the apparent 
impossibility of the ideal Toffoli and Fredkin unitaries un¬ 
der our desiderata (a)-(d). Although AC pulses were used 
in the proposed implementations with trapped ions and 
Bi donors, and for building circuits, these were merely a 
means to implement the time-independent TTfred and T^tof 
in appropriate rotating frames. This only results in extra 
relative phases between the computational basis states that 
do not matter for classical computing. White our gates 
minimize the dissipation in computation, the study of dis¬ 
sipation white pulsing and measuring and the possibility of 
room temperature realizations lls^ is kept for the future. 
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Interaction picture and the rotating wave approximation 

Consider a Hamiltonian that contains some fast oscillat¬ 
ing terms of the form V cos cot, where V is time indepen¬ 
dent and Hermitian. If a; > ||H||, then the evolution opera¬ 
tor can he approximated hy (see e.g. llh^ ) 


= l + (12) 

To represent this, we will use the notation iT (t) = Hq (t) + 
O ■ We now derive two identities that are useful in 

the main text. For this first identity, consider a Hamiltonian 
term of the form cos(a;a;t)cr^e*‘^‘^^*/^. Setting 


(jjx = ut gives 

= cos(u;f)c7^ = 

= -h -f 




c / ^ - - 

= — (cr^ -he a + e 


+a~ 


C 




c 


(13) 


provided that C < 2a;. For the second identity, consider 
a Hamiltonian of the form _|_ 

alal + erf +‘^='^21/2 

= + a^a^) + Jafa^ 

= + Jerfaf 

4J \ 

(14) 


= Jcrfcr| -F O 

\a;i - 012 
provided 4J < (a;i — a; 2 ). 


Thresholds for fault-tolerant classical computation 


We will estimate the error threshold required to imple¬ 
ment a fault-tolerant classical Toffoli gate, using the sim¬ 
plest classical code, the 3-hit repetition code. We will use a 
similar analysis as in ifT^Ihlll . which is not a very rigourous 
analysis hut will give a rough idea of the kind of classical 
threshold we will need to achieve with the Toffoli gate. A 
simple classically fault-tolerant Toffoli circuit can he con¬ 
structed as follows: 



Syndrome Recovery 


Where the syndrome -i- recovery step is 
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Since the code is robust against single classical (bit-flip) er¬ 
rors, a rough expression for the fault tolerant gate threshold 
can be found by considering the ways in which 2 or more 
errors can occur. Two errors can be output by a single fault- 
tolerant Toffoli gate in the following ways: 

1 There are incoming errors in two or more of the in¬ 
puts, which can propagate through the gate to create 
two errors on one encoded qubit. Incoming errors 
would be from a single error exiting from a previ¬ 
ous syndrome / recovery step, which could happen at 
any of six points in the syndrome / recovery process. 
There are ^ 6*2 = 3 ways that two errors could be in¬ 
coming, so if the error of a bit flip in the gate is p the 
total probability is 3 x ( 6 p)^ = 108p^. 

2 There is an incoming error in one of the inputs, and 
an error in the first round of Toffoli gates. Incoming 
errors have probability 6 p, and the probability of one 
error in any of the three gates is 3p. There are 3 ways 
this can happen so the overall probability is 3 x 3p x 
6 p = 56p^. 

3 Two errors occur during the Toffoli gates. There are 
^C 2 = 3 ways this can happen, so the total probabil¬ 
ity is 3p^. 

4 One failure occurs during the Toffoli gates (three 
ways this can happen), and one during one of the syn¬ 
drome gates (3 ways this can happen). Total proba¬ 
bility for all three encoded qubits is 3 x 3p x 3p = 
27p2. 

5 Two failures occur during one of the syndrome gates 
(^C 2 ways this can happen per encoded qubit). Total 
probability is 3 x 3p^ = 9p^. 

6 One failure occurs during the syndrome gate and one 
during the recovery. Total probability 3 x 3p x 3p = 
27p2. 

7 Two or more failures during recovery; total probabil¬ 
ity 9p^. 

In total, this gives a rough probability of (108+56+3+27+ 
9 + 27 + 9)p^ = 239p^ of two errors occurring undetected 


in the fault-tolerant Toffoli gate. Thus following concate¬ 
nation of this encoding (see e.g. ifT^ ) we would expect er¬ 
rors in the gate below the threshold p < ^ 5 x 10“^ 

to give a fault tolerant Toffoli gate. 

Parameter set for Fredkin Gate with Bismuth donors 

To satisfy the conditions of the Fredkin gate, J 12 = 
2(7l/| - 7l/|) and JUj^ “ 1 ) = W “ 

For the minimal gate time (n = 1), this would give 
J 23 = -^712, or alternatively J23 could be tuned to dif¬ 
ferent fractions of J 12 by altering m and n. For exam¬ 
ple, we can achieve J 23 = Ji 2 (l + 10 “®) by setting 
n = 675 and m = 2340. The resulting gate time would 
be 675 X 0.5ns ^ 0.5ps = 10 “®T 2 , so such tuning would 
still not lead to large decoherence errors. 



FIG. 4. Adding a fourth qubit to the original setup, in order to 
relax the constraints on Jzz and provide a control to turn the gate 
on and off. 

With J 12 = — AI 2 ), the error term in is 

0 ( 2(/2 — 7|)/(7i — ff))’ which is minimal when /f = 
— /2 = /3 = To decrease this error, a fourth qubit 

could be included, as shown in Fig. |4l Adding this qubit 
E with Ising coupling Je 2 effectively adds another local 
magnetic held to qubit 2 and so changes the resonance con¬ 
dition to J 12 + Je 2 = 0 J 3 — L 02 (assuming that this qubit 
is set in the |0) state). Note however that this coupling also 
has an error associated to it, and so the optimal situation 
is with Je 2 = J 12 = {AI^ — AI 2 ) = 1.475GHz. This 
results in an error of around 11%. Adding yet another con¬ 
trol qubit E' with coupling Je '2 and choosing couplings 
such that J 12 = Je 2 = Je '2 = ^{AI§ — AI 2 ) results in 
errors of around 7% (adding any more becomes unrealistic 
as the control qubits may begin to interact significantly). 
The gate time in this situation would be ^ 1.5ns which is 
still significantly smaller than T 2 . 

Additionally, setting J 12 = {ui 2 — W 3 ) is not straightfor¬ 
ward, since the donor position is not continuously tunable. 
There are also oscillations in the exchange interaction that 
depend on the separation between donors and the orienta¬ 
tion of the donors relative to the crystal llh^l . although these 
oscillations can be minimised if donors are aligned along 
the [100] axis, or if strain is applied 1^ . Atomically- 
precise positioning of the donors is possible 1 ^ 1 ^ . 
so a possible approach is to position the donors at separa¬ 
tions of around 15-20nm such that J 12 ~ (a ;2 — ^ 3 ) and 
then use a magnetic held gradient to tune (jJ 2 — UJ 3 closer 
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to Ji 2 (this magnetic field would also decrease the error in 
H^]). Since it would he possible to bring a magnetic tip 
close to the sample, large ma gnet ic field gradients of up to 
10^ Tm“^ would be possible 11^ so that a ;2 — ws could be 
tuned by up to ±0.5 GHz. Alternatively, we could adopt 
the method in and use electric gates to tune the inter¬ 
donor couplings and increase hyperfine interactions, how¬ 
ever this might introduce extra noise due to charge fluctua¬ 
tions. 

We could go further and use an extra qubit as an on/off 
switch for the gate, which could be useful if we wish to 
concatenate several of these gates together. Consider the 
setup in Fig. |4] that has one additional external qubit, la¬ 


belled qubit E. For qubit E to act as a control, we just need 
to make sure that the resonances of qubits 2 and 3 only 
match when qubits 1 and E are in the |1) state, and are 
very different otherwise so that the Heisenberg coupling 
becomes effectively an Ising coupling. By finding the res¬ 
onance energy of qubit 2 under different settings of qubits 1 
and E, and choosing J 12 = J 2 E, we find that the following 
conditions must be satisfied 

A + 2 Ji2 ^ J 23 ; A ^ J23, A — 2 Ji2 — 0 (15) 

Altogether these mean that the conditions under which this 
on/off switch would work is A ~ J 12 , and A, J 12 ^ J 23 . 


